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ABSTRACT 



Context. The physics of the pulsar inner magnetosphere remains poorly constrained by observa- 



tions. Although about 2000 pulsars have been discovered to date, only little is known about their 
emission mechanism. Large vacuum gaps exist in the magnetosphere and a non-neutral plasma 
fills partially the neutron star surroundings to form an electrosphere. 

Aims. In a previous work, we showed that the differentially rotating equatorial disk in the pulsar's 
electrosphere is diocotron unstable in the non-relativistic regime. In this paper, we extend these 
results and study the relativistic and electromagnetic stabilisation effects by including the mag- 
| netic field perturbation and allow for relativistic speeds of the guiding centre, in a self-consistent 

manner. We use the electric drift approximation, valid for low-density plasmas. 
Methods. We linearise the coupled relativistic cold-fluid and Maxwell equations in the elec- 
tric drift approximation. The non-linear eigenvalue problem for the perturbed azimuthal electric 
field is solved numerically with standard technics for boundary value problems like the shooting 
method. The spectrum of the relativistic diocotron instability in a non-neutral plasma column 
, confined between two cylindrically conducting walls is computed. 

Results. For low-speed motions, we recover the eigenfunctions and eigenspectra of the non- 
relativistic diocotron instability. Our algorithm is also checked in the relativistic planar diode 
geometry for which an analytical expression of the dispersion relation is known. As expected, 
when the relativistic and electromagnetic effects become significant, the diocotron instability 
tends to stabilise. In cylindrical geometry, for some special rotation profile, all azimuthal modes / 
are completely suppressed for sufficiently relativistic flows. However, for the profile relevant to 
the electrosphere, depending on the exact rotation curves, the growth rates can either significantly 
decrease till they vanish or persist for moderate /. 

Conclusions. The non-neutral plasma flowing in the pulsar electrosphere approaches the speed 
of light when reaching the light-cylinder. Therefore, relativistic and electromagnetic effects are 
important. They are capable to completely suppress the diocotron instability. Nevertheless, results 
are sensitive to the tail of the rotation curves. Therefore, particle diffusion across the magnetic 
field due to the diocotron instability only works efficiently close to the neutron star surface. 



Key words. Instabilities - Plasmas - Magnetohydrodynamics (MHD) - Methods: analytical - 
Methods: numerical - Pulsars: general 
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1. INTRODUCTION 



The detailed structure of charge distribution and electric-current circulation in the closed mag- 
netosphere of a pulsar remains poorly understood. Although it is often assumed that the plasma 
fills the space entirely and corotates with the neutron star, it is on the contrary very likely that it 
only partly fills it, leaving large vacuum gaps between plas ma-filled regions. The existence of suc h 



gaps in aligned rotators has been very clearly established by 



Krause-Polstorff & Michel 



( 1985a b). 



Since then, a number of differ ent numerical a pproaches to the problem have confirmed their conclu 
sions, including some work by Rvlov (1989), 



Thielh eim & Wolfstellerl dl994l) . 



Shibata(1989) 



Zachariadesl dl993), 



Neukirch(1993), 

Spitkovskv & Aronsl d2002l) . and by ourselves dPetri et alj2002b ). 



This conclusion about the existence of vacuum gaps has been reached from a self-consistent so- 
lution of Maxwell's equations in the case of the aligned rotator. Moreover, 



Smith et al 



(2001) 



have shown by numerical modelling that an initially filled magnetosphere like the Goldreich- Julian 
model evolves by opening up large gaps and stabil ises to the pa rtially filled and pa rtially void so- 



lution found bv lKrause-Polstorff & Michell dl985al) . and also by 



Petri et al 



(2002b). The status of 



models of the pulsar magnetospheres, or electrospheres, has recently been critically reviewed by 



Michell ( 12005b . A solution with vacuum gaps has the peculiar property that those parts of the mag- 



netosphere that are separated from the star's surface by a vacuum region are not corotating and so 
suffer differential rotation. 

This raises the question of the stability of such a charged plasma flow. The differential rotation 
in the equatorial, non neutral disk indu ces the so-called diocotron and magnetron ins tabilities that 



are well known to plasma physicists dOneil 



1980; 



Davidson 



1990; 



O'Neil & Smith 



1992). In the 



inner parts of the magnetosphere, far from the light cylinder, the instability reduces to its electro- 
static form, the diocotron instability. The linear d evelopment of the d ioco tron instabil ity of a thin 



differentially rotating charged disk was studied by 



Petri et al 



(2002a|) and 



Petril d2007) and shown 



to proceed at a growth rate c omparable t o the star's rotation rate. The non linear development of 
this instability was studied by 



Petri et al 



(2003), in the framework of an infinitely thin disk model. 
They have shown that the instability causes a cross-f ield transport of these charg es in the equatorial 
disk, evolving into a net out-flowing flux of charges. ISpitkovskv & Aronsl d2002l) have numerically 
studied the problem, and concluded that this charge transport tends to fill the gaps with plasma. 
The appeara nce of a cross-field elect ric current as a result of the diocotron instability has been 
observed bv lPasquini & Faiansl d2002l) in laboratory experiments in which charged particles were 
continuously injected in the plasma column trapped in a Malmberg-Penning configuration. 



Petri 



(2007) on the diocotron insta- 



The aim of this work is to extend the previous work done by 1 
bility by including electromagnetic and relativistic effects. Indeed, when the plasma approaches the 
light-cylinder, the guiding centre motion becomes relativistic and magnetic perturbations become 
significant. 

The relativistic aspect of the diocotron instab ility have already been investigated in the planar 



Davidson et al 



(11987 



1988). They clearly demonstrated the stabilisation due 



diode geometry by 

to electromagneti c effects. Applications to the planar magnetron geometry was investigated by 
Avres et al.l Jl992h . 



Send offprint requests to: J. Petri 
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In this paper we present a numerical analysis of the linear growth rate of the relativistic dio- 
cotron instability for a plasma column. The paper is organised as follows. In Sect. [2] we describe 
the initial setup of the plasma column consisting of an axially symmetric equilibrium between 
two conducting walls. We give several equilibrium profiles useful for the study of the diocotron 
instability in different configurations. In Sect. [3] the non-linear eigenvalue problem satisfied by 
the perturbed azimuthal electric field component is derived. The algorithm to solve the eigenvalue 
problem is checked against known analytical results in the non-relativistic cylindrical limit as well 
as in the relativistic planar diode geometry, Sect. [4] Then, applications to some typical equilibrium 
configuration are shown in Sect. [5] First we consider a plasma column with constant diocotron 
frequency. Next, we study the effect of the cylindrical geometry (curvature of the flow) and the 
transition to the planar diode limit. Finally, the stabilisation of the diocotron instability in a pulsar 
electrosphere is demonstrated when the flow becomes ultra-relativistic. The conclusions and the 
possible generalisation are presented in Sect. [6] 



2. INITIAL SETUP 



We study the motion of a non neutral plasma column of infinite axial extend along the z-axis, in 
the electric drift approximation. This approximation is well-suited for low-density plasmas satis- 
fying o»p <s a>l where a> p is th e plas ma frequency and a> c the cyclotron frequency. The geometric 
configuration is the same as in iPetril d2007l) . However, perturbations in magnetic field induced by 
the relativistic motion of the charges are now taken into account in a fully self-consistent manner. 
Therefore, we have to solve the full set of Maxwell equations. 

In this section, we briefly summarise the equilibrium conditions imposed on the plasma and 
give some typical examples of equilibrium configurations for specified velocity, density and electric 
field profiles. 



2.1. Equilibrium of the plasma column 

We consider a single-species non-neutral plasma consisting of particles with mass m e and charge q 
trapped between two cylindrically conducting walls located at r = W\ and r = W2 > W\. The 
plasma column itself is confined between R\ > W\ and R2 < W%- This allows us to take into account 
vacuum regions between the plasma and the conducting walls. We adopt cylindrical coordinates 
denoted by (r, <p, z) and the corresponding basis vectors (e r , e v , e z ). 

In the equilibrium configuration, the particle number density is n e (r) and the associated charge 
density is p e (r) = qn e (r). Particles evolve in a cross electric and magnetic field such that the 
equilibrium magnetic field is directed along the z-axis whereas the equilibrium electric field is 
directed along the r-axis. The electric field induced by the plasma itself, E p , is 

E p = E t e t . (1) 

The magnetic field is made of two parts, the first is an imposed external applied field, Bo, assumed 
to be uniform in the region outside the plasma column, and the second a plasma induced field, B p 



B = Bp + Bq = B z e z . 



(2) 
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Therefore, for azimuthally symmetric equilibria, the steady-state Maxwell-Gauss and Maxwell- 
Ampere equations satisfy 

-^irE t) = P -l (3 ) 
r or £o 

— = -HQPeV^. (4) 

In the electric drift approximation, particle trajectories are described by the guiding centre motion 
so that the flow velocity is only azimuthal and expressed as 

v = = -— e<p. (5) 

In order to avoid unphysical superluminal motion of the guiding centre, v < c, this model can only 
be applied for E r < c B z . In the opposite case, i.e. when E r > c B z , inertia effects become important 
and should be included in the equation of motion. The most general instability in a non-neutral 
plasma taking into account inertia and relativistic effects is the so-called magnetron instability. The 
full self-consistent response of the plasma to this instability will be studied in a forthcoming paper 
(removing the electric drift approximation, replaced by the relativistic momentum equation). In the 
case of pulsar magnetospheres, let us give an estimate of the distance to which this approximation 
remains valid. The drift motion, Eq. (0) applies roughly when the Larmor radius of the particles R g 
is less than the radius of their orbit R. Assuming relativistic speeds, with a guiding centre motion 
at the star rotation speed Q„, we get 

R Z^ = <R 

M B \q\B Vl -R 2 title 2 
Let R* and B* be the radius of the neutron star and the surface magnetic field respectively. For a 
dipolar vacuum magnetic field, we have 
R 3 

B(R) = B,^ (7) 

Condition Eq. (O can be cast into 

R 2 \q\B,Rl , 

= Y (8) 



yj\-R 2 & 2 Jc 2 m ec 
We introduced the parameter 



y=J^i (9, 

\ m e c 

Thus, we find that a rough estimate of the largest radius where the drift approximation remains 
valid is 



Y 2 R 4 L 

~ — A A 1 +4-4 - 1 (10) 

Applying to pulsar with period of the order one second, P — 1 s and magnetic fields of the or- 
der B, = 10 8 T, assuming a typical neutron star radius of R* = 10 km, we get 

fl L =4771/?* ; T = 5076fl L »/? L (11) 

Thus R mdx < Rl- The approximation falls down only very close to the light cylinder R^. For mil- 
lisecond pulsars, P = 1 ms and B„ = 10 5 T, the same conclusion applies because 

R L = 4.771 R* ; Y — 160531 R^ s> R^ (12) 
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We conclude that the electric drift approximation remains valid to good accuracy for electrospheric 
plasmas confined within the light cylinder. 

For a constant density profile in the plasma column with W\ = 0, the electric drift approxi- 
mation, Eq. ©, co rrespo nds to a circular motion at the diocotron frequency defined by, (see for 



instance 



Davidson 



(1990)), 



oj d = - — (13) 
2 u) c 

The relativistic plasma and cyclotron frequencies are given respectively by 

4 = (14) 

ym e e 



qB z 

ym e 



(15) 



where y = 1/ ^1 - v 2 /c 2 corresponds to the bulk Lorentz factor of the flow. 

We assume that the electric field induced by the plasma vanishes at the inner wall, at r = W\, 

i.e. 

EpiWO = 0. (16) 
Integrating Eq. ([3]) therefore gives for the electric field generated by the plasma, 

E P (r)=— f p e (r')r'dr'e T (17) 
sor J Wl 

For the magnetic field induced by the plasma, we solve Eq. dU with the boundary condition 
B z (Ri) - Bq. This simply states that the total magnetic field outside the plasma column has to 
match the magnetic field imposed by an external device. 

Any equilibrium state is completely determined by the following four quantities, the total radial 
electric field, E T , the total axial magnetic field, B z , the charge density, p e , and the azimuthal speed 
of the guiding centre, v v . Prescribing one of these profiles, the remaining three are found self- 
consistently by solving the set of Eqs. (|3), (|4]i and We give some typical examples in the next 
sections in which the velocity profile, the density profile or the electric field is imposed. 



2.2. Specified velocity profile 

Let us first assume that the velocity profile = r £2 is prescribed. This case is well-suited for 
the study of the pulsar's electrosphere in which the plasma is in differential rotation. Although 
the magnetic field is probably nearly dipolar in a pulsar, we assume a uniform mag netic field i n 
vacuum in order to remain fully self-consistent. We know from the previous analysis, ( Petrj|2007 ). 
that the growth rates are not very sensitive to the magnetic field structure. As already noticed in 
this work, the differential rotation is essential to the presence of the diocotron instability. The other 
equilibrium quantities, (E r ,B z ,p e ), are easily derived from £1 Indeed, inserting p e from Eq. © 
and E r from Eq. (O into Maxwell-Ampere equation ©, the magnetic field satisfies a first order 
ordinary linear differential equation 

^.tlB.L.n.t^U^ay (18) 

or r or c z or v ' 
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The Lorentz factor of the flow is 
1 

y 



rQ. 

P = — 

c 



(19) 
(20) 



From Poisson equation, Eq. ||3}, and the definition of the plasma frequencies, Eqs. (fT4b and ( TT3b . 
the charge density is found by 

> 2 „2 



Pe 



Til 

r dr 



(r 2 Q) = -y 2 



2D + r 



an 
~b7 



(21) 



The electric field is recovered from Poisson equation, Eq. or, according to the electric drift 
approximation, by Eq. (0 leading to the charge density 



Pe = so — tt (r i2 + r Q — — 

r Or y 1 or 

consistent with Eqs. (TT~8T > and ( f2TT >. 



(22) 



2.3. Specified density profile 

A simple and useful charge density profile to study the relativistic effect on the diocotron instability 
is given by a constant diocotron frequency in the whole plasma column, Eq. (fOT l, such that 

, Wi < r < R] 

oj-oir) = - co D = const , Ri < r < R 2 (23) 
, R 2 < r < W 2 

Using Eq. (I2TI 1. the rotation profile is deduced by integrating an ordinary differential equation for 
Q 

r 2 Q 2 ' 



Q + (L) D 1 



(24) 



with the boundary condition that Q vanishes at the inner plasma edge, O(/?0 = 0, to insure consis- 
tency with Eq. ( fToT ). Knowing £2, the same procedure as in the previous subsection for a specified 
velocity profile is applied, i.e. the magnetic field is calculated according to Eq. (1 1 8t> . 



2.4. Specified electric field 

It is also possible to specify the equilibrium radial electric field. An interesting case is given by 



EJr) 







-fin 



sinho'(r - Ri) 



, Wi < r < R y 
, Ri<r<R 2 



cosha(R 2 - R\) 

R 2 tanha(R 2 ~Rd 
-B — — , R 2 < r < W 2 



(25) 



a is a constant useful to adjust the maximal speed of the column at R 2 
v max = tanha (R 2 -R l ). 



(26) 



The equilibrium electric profile, Eq. (|25l ), enables us to investigate the influence of the cylindrical 
geometry compared to the planar diode geometry. Indeed, in the limit of small curvature of the 
column, i.e. when R 2 - R\ <K R\, the eigenvalue problem in cylindrical geometry reduces to the 
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planar diode case. The charge density is found by Eq. (0 whereas the magnetic field is solution of 
an ordinary differential equation 

dr r 2 c 2 dr 

Finally, the velocity is given by the electric drift approximation, Eq. (0. 



dB 2 1 d , , 9 x 

1^ = 7T-2TM E <) ^ 



3. LINEAR ANALYSIS 

In this section, we show how to derive the eigenvalue problem for the diocotron instability in the 
relativistic regime. Let us start with the full set of non-linear equations governing the plasma. 

3.1. Equations of motion 

Each particle evolves in the self-consistent electromagnetic field partly imposed by an external 
device and partly induced by the plasma itself. The motion of the plasma column is governed by 
the conservation of charge, the electric drift approximation, and the full set of Maxwell equations, 
respectively, 

^+div(p e v) = (28) 
dt 



E AB 



(29) 



B 2 
dB 

rotE = -— (30) 
dE 

rotB = fioj + so//o — (31) 
dt 

divE = — (32) 
so 

divB = (33) 

For a non-neutral plasma, the current density is related to the charge density by 

j=p e v (34) 

We apply the standard linear perturbation theory. All perturbations of physical quantities X like 
electric field, density, and velocity components, are expressed by the expansion 

X(r,ip,t) = X(r)e i< - l,p - l0t) (35) 

where I is the azimuthal mode and co the eigenfrequency. Therefore, we do not take into account 
the finite thickness of the disk. The whole cylinder moves as one block in the vertical direction. 
However, because the disk of the pulsar electrosphere has a finite vertical dimension, denoted by H, 
these modes should also depend on the vertical coordinate z. It is possible to extend this cylindrical 
model to modes depending on z by introducing the following expansion 

X(r, ip, z, t) = X(r) e i(kz+l (36) 
The technic is similar to the one used to investigate wave propagation in accretion d isks of finite 



thickness. Several methods could be applied. For instance. iNowak & Wagoned dl 99 II) used a series 
expansion in z/H of the Lagrangian displacement of a fluid element in order to study trapped 
adiabatic oscillations in accretion disks. Looking for solutions that are separable in the radial and 
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vertical direct i on is another mean to seek for the perturbations. This was done for instance by 



Okazaki et al 



Katold 19891) for 



( 1987) for trapped oscillations in relativistic accretion disks and by 
low-frequency corrugation w aves in an i s otherm al relativistic accretion disk. To get insight in the 
properties of the oscillations. iPerez et al.l d 1997b performed a radial and vertical WKB analysis. 

In our case, for a small vertical extension of the disk, k H » 1 , a series expansion in z/H would 
be appropriate. Indeed, regions with z » H will not contribute significantly to the electromagnetic 
field because of the oscillating term e' kz . However, this more general study is left for future work. 



3.2. Linearisation 



We study the stability of the plasma column around the equilibrium mentioned in the previous 
section. An expansion to first order for the electromagnetic field around the equilibrium (E°,B Z ) 
leads us to 



E = E° + 6E 

B z = Bl + 8B Z 

and the same for the charge and current density 
j = f + Sj 

Pe = Pe + ^Pe 

Linearising the set of Maxwell equations, we have 

1 d I 6p e 
— —{ro E T ) + i - SE^ = 

r s 



r dr 
1 d 



1 



-{rdE^) - i - 5E r 
r or r 



i to 5B 7 



I o> 

i-6B z = no 6j t - i -r 8E t 

r c L 

d 0) 

-—5B Z = fi 6j v - i -r 6E V 

or c 1 

The current density perturbation is 

6j t = p e 8v T 

Sj<p = ^Pe Vip + Pe 6V V 



(37) 
(38) 

(39) 
(40) 

(41) 
(42) 
(43) 
(44) 

(45) 
(46) 



It is convenient to introduce a new function (f> related to the azimuthal electric field by dE^ = -il<f>/r 
(we emphasise that this function is not the scalar potential from which the electric field could be 
derived from, it is just a convenient auxiliary variable) such that the electric and magnetic field 
become 



5E r = -K(r,u>) 



— - i no a) -r p e 6v T 

or l l 



tor 

6B Z = — K(r, oj) 



1 d(f> ,po 
-j^--i—p e 6v r 

c z or oj 



We introduced the function 
1 



K(r, oj) — - 

1- (cor lie) 2 

Maxwell-Gauss equation, d32l . is therefore written 



1 d 



r dr 



dr 



r- 



u> d l r 



- -r- \rK{r, cS) — - <f> = + ip Q — — -r <r, u>)p5v r 



£() 



r dr\l 2 



(47) 
(48) 

(49) 
(50) 
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From the continuity equation, d28l l, we get 
1 



1 d I 
- T-0"Pe <5v r ) + i-p6v (p 
r or r 



i{to-lCl) 
The electric drift approximation gives 



6v r = 
5v w = 



6E a 



I 



B z r B z 

E r 8B Z - 6E r B z 



(51) 

(52) 
(53) 



After some algebra, the eigenvalue problem for the relativistic diocotron instability in cylindrical 
geometry is written 



I d_ 

r dr 

with 



rn{r,io) — 
or 



I 2 

-r[l +x(r,a>)], 

r L 



/< 



r(to-l£l) 



K(r, to) 



D,r tor 
c Ic 



d_ 

dr 



tLl p 



U V I r V 

X(r, to) = K(r, to) — — 

CO c \lcl 



ftol 



h 2 K(r, to) — 

to c l 



(54) 



(55) 



The eigenvalue equation ( l54l l is very general. It describes the motion of small electromagnetic 
perturbations around the given equilibrium state, Eqs. (01 and in the electric drift approximation 
Eq. (0. Many aspect of the relativistic diocotron instability can be investigated with this eigenvalue 
equation. In order to solve the eigenvalue problem, boundary conditions need to be imposed at the 
plasma/vacuum interface. They play a decisive role in the presence or absence of the instability. 
How to treat these transitions between plasma and vacuum is discussed in the next subsection. 
N ote th at in the non-relativistic limit, the generalised linear eigenvalue problem, Eq. (20) in 



Petril d2007l) . is recovered 



\_d_ 
r dr 



dth 
Tr 



/< 



d 



r(to-lCi) dr 



I . ,2 



to c 



(56) 



This is the standard eigenvalue problem for the non-relativistic diocotron instability in cylindrical 
geometry. 



3.3. Boundary conditions 

In laboratory experiments, the plasma is usually confined between an inner and an outer conducting 
wall. However, in pulsar electrospheres, no such outer device exists to constraint the electric field 
at the outer boundary. Radiation from the plasma could propagate into vacuum to infinity, carrying 
energy away from the plasma by Poynting flux. To allow for this electromagnetic wave production 
by the instabilities studied in this work, the outer wall is removed. The electromagnetic field is 
solved analytically in vacuum and matched to the solution in the plasma at the plasma/vacuum 
interface located at r = R2- First we discuss the situation in which an outer wall exists and next 
consider outgoing waves. 



3.3.1. Outer wall 



When vacuum regions exist between the plasma column and the walls, special care is required at 
the sharp plasma/vacuum interfaces. Indeed, the right-hand side of Eq. d54l ) then involves Dirac 
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>\l > ULldlDL ILIllHUni (.'"'' 

words, its derivative is 



distribution functions 6(r) because the function toz/u) c is discontinuous at R\ and R 2 . In other 



d_ 
~dr 



5(r-R x ) - 
or 



( , ,2\ 



r=R, 



d 


i , l2^ 


d 


( , ,2\ 


dr 


LOr 


+ d~r 


W c 1 




\ L I 


r=R 2 





(57) 



regular 



where || re guiar means the regular (or continuous) part of the derivative, i.e. which does not involve 
distribution functions. It vanishes in the vacuum regions, r < R\ and r > R 2 . Therefore, the first or- 
der derivative of <p is not continuous at these interfaces. To overcome this difficulty, we decompose 
the space between the two walls into three distinct regions: 

- region I: vacuum space between inner wall and inner boundary of the plasma column, with the 
solution for the function (p denoted by (pi, defined for W\ < r < R\ ; 

- region II: the plasma column itself located between R\ and R2, solution denoted by (pn, defined 
for Ri < r < R 2 ; 

- region III: vacuum space between the outer boundary of the plasma column and the outer wall, 
solution denoted by (pm, defined for/?2 < r <W 2 - 

In regions I and III, the vacuum solutions should satisfy the required boundary conditions, 4>i(W\ ) = 
and (p m (W 2 ) = 0. 

The jumps in the derivative d(p/dr at each interface are easily founded from Eq. d5Tb . At R\, the 
jump is 

Q(7?i)/?i ojRA o$(Ri) 



dr dr R x (oj-IQ.(Ri)) 



Ic J u) Q (R\) 



Similarly, at the outer interface at R 2 , we obtain, 



90m , <30 u 
-^—(Ri) — H - ("2) 
or or 



1<P(R 2 ) 



R 2 (cj-IQ(R 2 )) 
3.3.2. Outgoing wave solution 



£l(R 2 )R 2 cjR 2 \ (^(Ri) 



lc j OJciRl) 



(58) 



(59) 



Because of the wall located at r — W 2 , the outer boundary condition (pm(W 2 ) = enforces 
E,f(W 2 ) = 0. It therefore prevents escaping waves from the system due to the vanishing outgo- 
ing Poynting flux, E ip B z /fio = 0. In pulsar magnetospheres, no such wall exists. So, in order to 
let the system produce outgoing electromagnetic waves, we remove the outer wall in this case and 
solve the vacuum wave equation for cp which then reads 

1 d 



dr 



ric{r,aj) — 
dr 



I 2 

~^cp = 



This equation can also be derived directly from the vector wave equation 

AE ~7 2 W = 



(60) 



(61) 



projected along the axis. To find the right outgoing wave boundary conditions, it is therefore 
necessary to solve the vector wave e quation in cylin dric al coordinates using vector cylindrical har- 



monics as described for instance in 



Stratton 



dl941l) and lMorse & Feshbachl dl953l) . The solutions 
for the function <p to be an outgoing wave in vacuum outside the plasma column and which vanishes 
at infinity is given by (region III with W 2 = +00) 
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11 



where the cylindrical outgoing wave function is given by Hi dStratton 



1941), the Hankel func- 
tion of first kind and of ord er I related to the Bessel functions by H[{x) - Ji(x) + i Y[{x), 



dAbramowitz & Stegunl 11965). The prime ' means derivative of the function evaluated at the 
point given in parentheses. K is a constant to be determined from the boundary condition at R 2 . 
Eliminating the constant K, we conclude that the boundary condition to impose on (p is 



„,(uRi\ ojR 2 „(uR 2 \] iW w *2\ <^>m 



(63) 



The boundary conditions expressed in region II for <pu are found by replacing 0' m (R 2 ) from Eq. 
and recalling that <p is continuous, therefore (f>m(R 2 ) = (puiRi) = QiRi)- We find 



a>R 2 
c 



<P(R 2 ) - 
1<P(R 2 ) 



(64) 



dr^ Rl) R 2 (oj - m(R 2 )) 



1 - 



Cl(R 2 )R 2 uR 2 \ ^1{R 2 ) 



Ic ) Ll> c (R 2 ) 



= 



3.4. Algorithm 

The eigenvalue problem, Eq. (l54l . is solved by standard numerical technics. We have implemented 
a shooting method as follows. 

First, we guess a good initial value for the eigenvalue co. Then, the ordinary differential equa- 
tion, Eq. (l54l . is integrated numerically from W\ to W 2 with a fifth-order Runge-Kutta or a Bulirsch- 
Stoer method. More precisely, at r — W\, the initial conditions are, (p = and d<p/dr - 1 . Integration 
is performed in region I until the first vacuum/plasma interface is reached at r = R\. There, the 
first order derivative in <p is subject to a discontinuity given by the jump in Eq. d58l l. Knowing 
d(j)ildr(R\), we deduce d<pu/dr(Ri). Therefore, the integration is continued in region II with the 
initial conditions, (/>n(Ri) = <pi{Ri) (because <p is continuous) and d<pu/dr(R 2 ) until the second 
vacuum/plasma interface is reached at R 2 . The first order derivative in <p is now subject to an- 
other discontinuity given by the jump in Eq. d59b . Knowing d(pn/dr(R 2 ), we deduce d<f>m/dr(R 2 ). 
Integration is continued in region III with the initial conditions, (pmiRi) = 0nC?2) an d d<f>m/dr(R 2 ) 
until the endpoint W 2 . At the end of the process, the function 0ni does not necessarily satisfy the 
desired boundary conditions. Indeed, the eigenvalue is found whenever the function <p at the outer 
wall vanishes <pm(W 2 ) = 0. Finding <±> is therefore equivalent to finding the root of </>m(W2) with 
respect to the eigenvalue a>. 

For the pulsar electrosphere, the situation is very similar except that no calculation is performed 
in region III. The boundary condition for outgoing waves is applied at R 2 , see Eq. d64i >. Actually 
for pulsars, we compare both boundary conditions. 



4. Algorithm check 

In order to check our algorithm in different configurations, we compute the eigenvalues for both 
a non-relativistic cylindrical plasma column and a relativistic planar diode geometry. For some 
special density profiles, the exact analytical dispersion relations are known and used for comparison 
with the numerical results. 
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4.1. Non-relativistic plasma column 

In cylindrical geomet ry, an exact analy tical solution for the dispersion relation can be found in the 



1990). We use these results to check our algorithm in cylindrical 



non-relativistic case, (iDavidson 
coordinates. 

The magnetic field is constant and uniform in the whole space, B z (r) = Bq. We do not need 
to solve Maxwell-Ampere equation because the magnetic perturbation is neglected in the non- 
relativistic limit. The particle number density and charge density are also uniform in the whole 
plasma column such that 



, Wi < r < Ri 

Pe(r) - \ po - const , R\ < r < R 2 
, R 2 < r < W 2 

Solving Maxwell-Gauss equation d32b , the equilibrium radial electric field is 



(65) 



EAr) = 





Por 
2 e 

P() R 2 



Rp 



Wi < r <Ri 
Ri < r < R 2 

R 2 < r<W 2 



(66) 



2 e r 

and the corresponding electric drift speed in the plasma 



2\ 



R 



(67) 



The diocotron frequency is constant, wd = pqI2eqBq = const. The solutions to the dispersion 
relation for this particular case are 



u> = -- 



Ll>D 



b,+ yjb^-4ci 



The coefficients b\ and c\ are given by 



b, = ll 



1 



1-1* 

Ri 

1-1* 

Ri 



Rj_ 

w 2 



w 2 ) 
w 2 ) 



i-i^ 



w 2 



21 



1 - 



(68) 



(69) 



1 



^7 



W 2 



(70) 



A sample of eigenvalues is shown in Tab. [T] for W\ — 1 and W 2 = 10 and different aspect ra- 
tios, d\ - R1/W2 and d 2 = R 2 /W 2 . The exact analytical solution, Eq. d68b . is compared with the 
numerical solution. The relative errors in the real and imaginary part of the eigenvalues are 
Re (oj) - Re (awO 



film 



Re (Wexact) 

Im (o>) - Im (w exact ) 



< 3.2 x 10" 



< 1.1 x 10" 



(71) 
(72) 



Im (w exact ) 

The precision is excellent, it reaches 10 digits. Our algorithm computes quickly and accurately the 
eigenvalues in cylindrical geometry with vacuum gaps between the plasma column and the walls. 
The eigenvalues obtained in this example are good initial guesses to study the relativistic problem 
in the low speed limit. Next, we turn to the relativistic planar case. 
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mode / 




d 2 


^num 








2 


0.4 


0.5 


-3.772e-01 + 7.176e-02 


i 


2.046e-ll 


3.148e-ll 


3 


0.4 


0.5 


-5.456e-01 + 2.267e-01 


i 


3.330e-15 


7.352e-15 


4 


0.4 


0.5 


-7.216e-01 + 2.988e-01 


i 


7.845e-12 


1.855e-12 


5 


0.7 


0.9 


-1.147e+00 + 5.787e-02 


i 


6.050e-12 


2.618e-ll 


7 


0.6 


0.7 


-9.315e-01 + 3.307e-01 


i 


1.006e-ll 


3.832e-12 



Table 1. Numerical eigenvalues w num and relative errors for the non-relativistic plasma column 
for different mode / and different aspect ratios, d\ = R1/W2, and d.2 - R2IW2 with W\ — 1 and 
W2 - 10. The precision of the computed eigenvalues reaches 10 digits. 



4.2. Relativistic planar diode geometry 



In the relativistic planar diode geometry, the dispersion relation is also known analytically in the 
long wavelength limit. Thus it is an interesting case to check our algorithm for the relativistic 
diocotron instability. For com pleteness, we recall the main results. For a detailed discussion, see 



Davidson et al 



(11987 . 



1988). The plasma is drifting in the y-direction at a speed V y and located 



between x — X\ and x — X2. The cathode is located at x = and the anode at x = d. We use 
Cartesian coordinates (x,y,z). The eigenvalue problem in Cartesian coordinates for the diocotron 
instability in a relativistic planar diode geometry is 



dx 2 



,2 \ 



1 



k 2 c 2 k 2 c 2 



Co - k V v 



Vy OJ \ diL>Y) { 



c kcl dx 



(73) 



where k is the wavenumber. In this paragraph, because of the Cartesian geometry, the diocotron 
frequency is defined as (no factor 1 /2) 

, >2 



OJ£>C - 

OJr 



The applied external magnetic field is constant and uniform outside the plasma layer, B z (r) 
for x > X2. The particle number density is chosen such that 



n(x) 
y(x) 





»0 



Wi < x < Xi 



— = const , X\ < x < X2 

70 

, X 2 < x < W 2 



(74) 
= B 

(75) 



The Lorentz factor is y — 1 / ^1 - Vy/c 2 . For this particular choice of density profile, the diocotron 
frequency is constant throughout the layer cross section, wd c = const. Moreover, the drift speed is 



V y - c tanh [wd c (x - Xi)/c] 



(76) 



In the long-wavelength perturbation limits, corresponding to k d <K 1 , the dispersion relation reads 

- 1 



Re(oj) 
k c 

lm(co) 
kc 



A; 

sinh 6» + 2 — 6»cosh6» 
Ab 




cosh 6 + 



e 



Ab sinh 



(A; cosh 6 + A ) 



sinh 2 9 
~9~ 2 ~ 



cosh 6 + 



Ai, sinh 9 



(A/ cosh 6 + A ) 



where the quantities are defined by 

A,- = X1/W2 



(77) 



(78) 



(79) 
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A h = (X 2 -X 1 )/W 2 (80) 

A = (W 2 -X 2 )/W 2 (81) 

p = 2 VaTAo/A,, (82) 

= oj Dc (X 2 -X 1 )/c. (83) 



Results for the eigenvalues are shown in Figs.Q]and[2]in the long wavelength limit and compared 
with the exact dispersion relations, Eqs. ( fTTT i and ((78). The relativistic diocotron instability for 
arbitrary wavelengt h is shown in F igs.[3]and[4] Our numerical results agree with very good accuracy 



to those found by 



Davidson et al 



(11987 . 



1988b who directly solved for the dispersion relation in 



the general case, i.e. for arbitrary kd. 




Fig. 1. Imaginary part, on the left panel, and real part, on the right panel, of the eigenfrequencies, 
(in black dots), compared to the analytical exact dispersion relation in the long wavelength limit 
Eq. f77) and dTBl . (red curves). The parameters are A = 0.3 and A, = 0.3, 0.4, 0.5, 0.6. 




Fig. 2. Imaginary part, on the left panel, and real part, on the right panel, of the eigenfrequencies, 
(in black dots), compared to the analytical exact dispersion relation in the long wavelength limit 
Eq. (|77]l and {78}, (red curves). The parameters are A, = 0.3 and A = 0.3, 0.4, 0.5, 0.6. 




Fig. 3. Imaginary part, on the left panel, and real part, on the right panel, of the eigenfrequencies 
for arbitrary wavelen gth, kd. The parameters are A,- = 0.3 and A ( > = 0.3. Check with Fig. 6 of 
Davidson et al.Nl988l) . 
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Fig. 4. Imaginary part, on the left panel, and real part, on the right panel, of the eigenfrequencies 
for arbitrary waveleng th, kd. The parameters are A,- = 0.3 and Aq = 0.6. Check with Fig. 7 of 



Davidson et al 



(1988). 



5. RESULTS 

We demonstrated that our numerical algorithm gives accurate results in the non-relativistic cylin- 
drical geometry as well as in the relativistic regime in Cartesian coordinates. In this section, we 
compute the eigenspectra of the relativistic diocotron instability in cylindrical coordinates, for var- 
ious equilibrium density, electric field, and velocity profiles. Application to pulsar's electrosphere 
is also discussed. 



5. 1 . Plasma column 

First, we consider the laboratory plasma in cylindrical geometry, confined by some external ex- 
perimental electromagnetic device. The external applied magnetic field and the density profile are 
specified as initial data. We use the simple profile presented in Sect. I2.3I 

In the non-relativistic limit, a very good first guess is given by Eq. d68l i. The influence of 
the relativistic and electromagnetic effects are then investigated by slowly increasing the maximal 
speed at the outer edge of the plasma column, v max = Ro Q{Ri)- 

In the non-relativistic flow, the growth rate of the diocotron instability only depends on the 
ratios d\ = R\IWi, dt = R2/W2 and w = W1/W2. In order to avoid variation of the growth rate 
or stabilisation of the instability purely due to geometrical effects, these ratios are kept constant 
while v max is increased. To this end, for instance, we increase W2 and adjust the other parameters 
R] , R2, W\ to fit the constant ratios conditions, w = est, d\ = est, d2 = est. 

Two cases are presented in Fig. [5] The first one has w = 0.1, d\ = 0.45 and di = 0.5, Fig.|5^), 
whereas the second one has w = 0.1, = 0.4 and ^2 = 0.5, Fig. \5&). The eigenvalues are 
conveniently normalised to the diocotron frequency con. For non-relativistic speeds, v max <c c, 
the eigenvalues of Sect. 14. H are recovered, see Table [2] The thinner the plasma layer, the larger the 
number of unstable modes, respectively 5 and 12 unstable modes. In both cases, the growth rate 
starts to be altered whenever v max /c > 0. 1 . We only labeled the first 5 unstable modes I = 2, 3, 4, 5, 6 
in order to avoid overloading the figure. 

In any case, for very high speeds, v max /c m 1, all the diocotron modes become stable because 
the growth rate vanishes, Fig. [5] The stabilisation process already observed in the relativistic planar 
diode is not altered by the cylindrical geometry. 
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Table 2. Eigenvalues of the plasma column for the density profile in Sect. 14.11 Comparison of the 
non-relativistic one w nre i and the relativistic one <x> re i in the low speed limit, v max <K c for w = 0.1, 
di = 0.4 and d 2 = 0.5. 



mode / Olnrel W re j 



2 


3.772986e-01 + 7.176435e-02 i 


3.773015e-01 +7.174847e-02 i 


3 


5.456744e-01 + 2.267241e-01 i 


5.456773e-01 + 2.267196e-01 i 


4 


7.216191e-01 +2.988911e-01 i 


7.216219e-01 + 2.988890e-01 i 


5 


9.004354e-01 + 3.118851e-01 i 


9.004382e-01 + 3.118852e-01 i 


6 


1.080114e+00 + 2.495296e-01 i 


1.0801 16e+00 + 2.495328e-01 i 




3 iTi 



Fig. 5. Stabilisation of the diocotron instability when relativistic and electromagnetic effects are 
included. The relativistic diocotron frequency is constant in the whole plasma column, a>£, = cste. 
The geometric aspect ratios are, w = 0.1, d\ = 0.45 and d% = 0.5 for Fig. a) and w = 0.1, 
d\ = 0.4 and di = 0.5 for Fig. b). Only the first unstable modes are labeled, I = 2, 3, 4, 5, 6 to avoid 
overloading the plots. Each coloured curve depicts a different azimuthal mode number I. 

5.2. Relativistic planar diode limit 

The influence of the curvature is also studied by taking the limit of the planar diode geometry. The 
curvature of the plasma column is then increased to investigate the evolution of the growth rates. 
The effective aspect ratio of the plasma layer is conveniently described by the parameter 

A = (84) 
R2-R1 

Using the equilibrium electric field profile indicated in Sect. 12.41 in the limit of small curvature cor- 
responding to large aspect ratio, A — > +00, the eigenvalue problem and equilibrium configuration 
is described by the relativistic planar diode. 

We show the evolution of the growth rate in the non relativistic limit, v max = 10~ 3 , Fig.|6^), and 
in the relativistic case, v max = 0.3, Fig. |6j)). The eigenvalues are conveniently normalised to the 
value of the diocotron frequency at the outer boundary of the plasma column, cjuiRi)- The aspect 
ratio has a drastic influence on the growth rate. For large values of A » 1, all unstable modes 
are stabilised, in both non-relativistic and relativistic flows. Note however, that in the relativistic 
regime, the mode / = 2 and I = 3 have already been stabilised, whatever the aspect ratio. 

When the aspect ratio A is increased, the fastest growing unstable mode is shifting to higher 
azimuthal numbers / whereas the lower azimuthal numbers / start to stabilise. Indeed, to compute 
the eigenvalues for different aspect ratios, we increase Ri and R2 while keeping the difference 
A = R2 - R\ constant. Assuming that the geometrical size of the perturbation fits into the layer 
thickness A, we get an estimate of the allowed azimuthal numbers such that 
2n 



(85) 



Therefore the fastest mode numbers are roughly I oc A and linearly growing with the aspect ratio. 
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Fig. 6. Effect of the cylindrical geometry on the growth rate of the relativistic diocotron instability. 
The growth rate are normalised to the diocotron frequency and plotted versus the aspect ratio A. 
Non-relativistic regime, v max = 1CT 3 in Fig. a) and relativistic speed, v max = 0.3 in Fig. b). Only the 
first unstable modes are labeled, I = 2, 3, 4, 5 to avoid overloading the plots. Each coloured curve 
depicts a different azimuthal mode number /. 

5.3. Electrosphere 

The electrospheric non-neutral plas ma, as a l ready proved in previous works by 



Krause-Polstorff & Michel (1985a) and 



Petri et al 



d2002bl) . is confined by the rotating mag- 



netised neutron star. The most important feature is the velocity profile in the plasma column. For 
simplicity, here, we assume that no vacuum gaps exist between the plasma and the walls, so that 
Wi = Ri and W2 = R%. The rotation profile is chosen to mimic the rotation curve obtained in the 
3D electrosphe re. To study the influence of the relativistic effects, we take the same profiles as 
those given in Petri (120071) . We remind that three different analytical expressions for the radial 



dependence of Q are chosen by mainly varying the gradient in differential shear as follows 
Q(r) = Q, (2 + tanh[a(r - r )] e^' A ) 



(86) 



The values used are listed in Table (f3j). The angular velocity starts from corotation with the star 



n 


a 


P 


ro 


SI] 


3.0 


5 x 10~ 5 


6.0 


n 2 


1.0 


5 x 10- 5 


6.0 


n 3 


0.3 


5 x 10~ 5 


10.0 



Table 3. Parameters for the three rotation profiles used to mimic the azimuthal velocity of the 
plasma in the electrospheric disk. 



D. = £2,, followed by a sharp increase around r — 6 for £2^2, and a less pronounced gradient around 
r — 10 for Sl j,. Finally the rotation rate asympto tes t wice the neutron st ar rotation speed, Fig. [7] 



The results of 



Krause-Polstorff & Michel 



Petri et al. 



(2002b) have been computed for 



d!985al) and] 

non-relativistic speeds. However, confinement of non-neutral plasmas by some external electro- 
magnetic fields is very ge neral and also applies when relativistic effects are taken into account, 



Tsang & Davidson! d 1 9860 - In case of relativistic motions within the electrosphere, we would ex- 



pect some quantitative changes in the rotation curves because of the relativistic factor appearing in 
the expression for the charge density (due to the current displacement) 
£1 B 



p--2s - 



(87) 



1 -r 2 Q 2 /c 2 

Nevertheless, differential rotation is still expected, even in this latter case. That is why we took 
typical rotation curves as those depicted in Fig. [7] a). We emphasise that the electromagnetic field 
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Fig. 7. Three choices of differential rotation curves in the plasma column for the cylindrical pulsar 
electrosphere, fii in red, Q.2 in green and Q3 in blue, Fig. a). The linear speed is depicted in Fig. b) 
in units of v max = R2 CliRi)- 



as well as the charge density resulting from the chosen rotation profiles are determined in a full 
self-consistent manner, as described in Sect. [2] Relativistic speeds are only reached in the extended 
part of the electrosphere, i.e. the outer part approaching the light cylinder. In Fig. 0b), the linear 
speed in the disk is plotted r £2(r)/v max , the highest speed v max < c corresponds to the largest radii r, 
i.e. the outer part of the disk. Nevertheless, in the region where strong gradient exists, around r — 6, 
the denominator of Eq. ( f87b remains close to unity. Indeed, in the differentially rotating part, we 
have r Q(r)/v max < 1/2 implying (1 - r 2 Q?/c 2 y l < 4/3, therefore the non-relativistic calculations 
are still valid in this part of the disk with an error less than roughly 20%. Actually, the exact rotation 
profile is not important in our study, we just want to demonstrate that the diocotron instability still 
exists in the relativistic regime. 



We generalise the study presented in Petri (2007) by including the relativistic effects, i.e. rela 



tivistic speed of the flow and electromagnetic field perturbation. 

We start with a non-relativistic rotation profile such that Q.* R2 <s c and slowly increase R2 (as 
well as W\, W2,R\ to maintain their ratio constant) in order to approach the speed of light for the 
maximal rotation rate of the plasma column. 

First, we computed the growth rates when an outer wall is presented. Second, we remove this 
wall, imposing outgoing wave solutions and finally compare both situations. 

5.3.1. Outer wall 

To remain fully self-consistent, we only consider an uniform applied external magnetic field. The 
maximum growth rates, normalised to the angular velocity of the neutron star, for the first rotation 
curve, Qi, for each mode /, is shown in Fig. [8^) for an extension from Ri = 1 to R2 — 20. Each 
coloured curve depicts a different azimuthal mode number /. However, they are not labeled because 
we want only to demonstrate the stabilisation of all modes. The precise value of the eigenvalues 
associated with a particular mode I therefore does not matter in our present study. The most impor- 
tant characteristic is the behaviour of the full set of unstable eigenvalues when the system becomes 
relativistic. 

The profile having the steepest gradient possesses the largest number of excited unstable modes 
because it corresponds to the case where the smallest scales appear, i.e. I large. In the non- 
relativistic limit, v max <k c, the l argest grow th rate, for / = 8 has a value of y max = 3.2. The 



spectrum is exactly the same as in Petri (2007). When the maximal speed is increased, the growth 



rates vary significantly but we do not observe the stabilisation effect except for the mode I = 1 
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which disappears for v max = 0.936. For relativistic speed, v max c, the other modes keep roughly 
their growth rate at their value for the non-relativistic instability. Several of them even increase. 

In order to demonstrate the tendency towards stabilisation in the electrosphere, we reduce the 
size of the plasma column. For instance, we choose the plasma extension from Ri — 1 to R2 = 10. 
The new growth rate are shown in Fig. [8})). The diocotron instability now tends to stabilise for all 
modes I. Close to the speed of light, the growth rates start to decrease significantly. Note that the 
I = 1,2 modes already disappeared, Fig.lSJ)). 

The second steepest profile possesses less unstable modes as we would expect due to the fact 
that only larger scale structures can emerge with this slope of the differential rotation, Fig. |9^). 
Here, the tendency to stabilise the diocotron instability is more evident. The mode / = 1 disappears 
as in the previous case. The other modes commence to show a significant decrease in their growth 
rate when approaching the ultra-relativistic limit. Here again, we checked that for low speeds, we 
recover the non-relativistic spectrum with good accuracy. For the narrower layer, the stabilisation 
is fully achieved for all the modes, Fig. [9j>). 

Finally, the third smooth profile has only four unstable modes, Fig. [TOh). The stabilisation 
effects are clearly seen for the modes I = 2 whereas the other modes are near to full stabilisation. 
Reducing the size of the electrosphere, here again we observe full disappearance of the diocotron 
instability, Fig.fTOb). 




9^ j£V t.t a.t □.■ 1 » a.i m.u □■■ 1 



Fig. 8. Evolution of the growth rates lm(a)) of the profiles £2i for increasing maximal speed of the 
column. The column is extending from Ri — 1 to R2 — 20 in Fig. a), and from Ri — 1 to R2 = 10 
in Fig. b). Each coloured curve depicts a different azimuthal mode number /, they are not labeled. 
The outer wall is located at 7?2- 



The proof of the relativistic stabilisation effect are therefore obvious for the profiles Q2 and O3. 
It is clearly seen that the diocotron instability is suppressed when the maximal speed approaches 
the speed of light with a steep decrease in growth rate for v max ^ c. These examples undoubtedly 
reveal the influence of relativistic and electromagnetic effects towards stabilisation of the non- 
neutral plasma in the pulsar electrosphere. 

To understand the change in the behaviour of the diocotron instability when simply changing 
the size of the plasma column, or equivalently, the dimension of the electrosphere, we need to 
understand how this instability is generated. The instability is launched from the resonance between 
the motion of the wave pattern related to the perturbation and particles in the flow. It is readily seen 
from the denominator of Eq. (1541 . The corotation radius r c satisfies 

Re{oS)=lQ.{r c ) (88) 

In order for the relativistic effects to play a significant role, the speed in the vicinity of the corotation 
region should be close to the speed of light. If the tail of the rotation curve is long, the outer edge 
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Fig. 9. Evolution of the growth rates lm(a)) of the profiles Q.2 for increasing maximal speed of the 
column. The column is extending from R\ — 1 to Ri — 20 in Fig. a), and from R\ — 1 to R2 — 10 
in Fig. b). The outer wall is located at Ri- 




Fig. 10. Evolution of the growth rates Im(o») of the profiles Q3 for increasing maximal speed of the 
column. The column is extending from R\ — 1 to R2 — 20 in Fig. a), and from Ri — 1 to Ro — 15 
in Fig. b). The outer wall is located at Ro- 
of the plasma attains speeds close to c while the speed near the corotation radius remains non- 
relativistic. Thus, the growth rates are not significantly affected by the relativistic effects. However, 
reducing the tail of the rotation curve allows the region near corotation to reach higher velocities. 
Therefore, the stabilisation starts to set in and the diocotron instability can be suppressed. 

5.3.2. Outgoing waves 

We performed a second set of calculations by removing the outer wall assumption and enforce 
outgoing electromagnetic waves propagation into vacuum. The boundary condition on the outer 
plasma/vacuum interface has been discussed in Sect. 13.31 

We use exactly the same configurations and rotation profiles presented in the previous section. 
The only change comes from the outer boundary condition, namely ipmi^i) = which is replaced 
by Eq. ( 164) . The results are shown in Fig. QT| [T2l and [T3] for the rotation curve Qi, Q2 and Q3 
respectively. Comparing both situations, the growth rates are not significantly affected by wave 
emission. Note however, that because energy is carried away by Poynting flux, the instability grows 
slowlier than in the previous case. 




Fig. 11. Same as Fig. <[8J but outgoing waves boundary conditions are applied. 
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Fig. 12. Same as Fig. (0 but outgoing waves boundary conditions are applied. 




Fig. 13. Same as Fig. ( TTOb but outgoing waves boundary conditions are applied. 



In Fig. [T4J we show a comparison between both boundary conditions for the profile Q2. 
It is clearly seen that the growth rates are relatively insensitive to the nature of the boundary. 
Nevertheless, there is a tendency to decrease the growth rates when outgoing waves are present. 




Fig. 14. Comparison of the outer wall and outgoing wave boundary conditions for the profile Q2, 
Fig 0b) and Fig[l2]b). The growth rates for outgoing waves are almost always slightly smaller. 



6. CONCLUSION 

We developed a numerical code to compute the eigenspectra and eigenfunctions of the diocotron 
instability including electromagnetic and relativistic effects. In the non-relativistic limit, we re- 
covered and confirme d the results obtained by another technic using pseudo-spectral method by 



Chebyshev expansion, Petri (2007). Unstable modes are computed for a uniform external applied 
magnetic field and arbitrary velocity, density and electric field profiles. Application to a plasma 
column as well as to the pulsar electrosphere have been shown. In both cases, the diocotron regime 
gives rise to instabilities with decreasing growth rates when the maximal speed of the flow becomes 
relativistic. Whereas the growth rates can be comparable to the rotation period of the neutron star 
in the non-relativistic limit, it is found that for special rotation profiles, the diocotron instability 
is completely suppressed in the relativistic regime. Including electromagnetic wave emission from 
the electrospheric plasma does not change drastically these conclusions. 

What therefore happens to the plasma in the vicinity of the light cylinder needs a more general 
treatment including inertia of the particles, because the plasma kinetic energy becomes comparable 
to the magnetic field energy density. At the light cylinder, the plasma will therefore be subject to the 
magnetron instability. It is the most general case, including relativistic flow, electromagnetic pertur- 
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bation and inertia of the particles. The study of the magnetron instability in a pulsar electrosphere 
is the aim of a forthcoming paper. 

Last but not least, the influence of finite temperature in the plasma on the diocotron or mag- 
netron instability would require a kinetic treatment of the stability via the Vlasov-Maxwell equa- 
tion. This is also left for future work. 
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